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Abstract 

Within a two-loop leading-log approximation, we review the prediction for the lightest 
Higgs mass (m/i) in the framework of constrained MSSM (CMSSM), derived from the 
naturalness requirement of minimal fine-tuning (A) of the electroweak scale, and dark 
matter consistency. As a result, the Higgs mass is predicted to be just above the LEP2 
bound, uih = 115.9 ± 2 GeV, corresponding to a minimal A ~ 17.8, value obtained from 
consistency with electroweak and WMAP {3a) constraints, but without the LEP2 bound. 
Due to quantum corrections (largely QCD ones for ruh above LEP2 bound), A grows « 
exponentially on either side of the above value of nth , which stresses the relevance of this 
prediction. A value rrih > 121 (126) GeV cannot be accommodated within the CMSSM 
unless one accepts a fine-tuning cost worse than A> 100 (1000), respectively. We review 
how the above prediction for riih and A changes under the addition of new physics beyond 
the MSSM Higgs sector, parametrized by effective operators of dimensions d=5 and d=6. 
For d=5 operators, one can obtain values nih < 130 GeV for A < 10. The size of the 
supersymmetric correction that each individual operator of d=6 brings to the value of 
mh for points with A < 100 (< 200), is found to be small, of few < 4 GeV (< 6 GeV) 
respectively, for M = 8 TeV where M is the scale of new physics. This value decreases 
(increases) by approximately 1 GeV for a 1 TeV increase (decrease) of the scale M. The 
relation of these results to the Atlas/CMS supersymmetry exclusion limits is presented 
together with their impact for the CMSSM regions of lowest fine-tuning. 
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1 Naturalness measures and predictions. 



After forty years of super symmetry, we are fortunate to currently witness its biggest real 
test at the Large Hadron Collider (LHC). It remains to be seen if the LHC experiments will 
prove the idea of low-energy (TeV-scale) supersymmetry, or just increase its scale to higher 
values. Negative searches for superpartners in the TeV region can ultimately question it as 
a natural solution to the hierarchy problem. This solution resides on quantum cancellations 
between partners and superpartners, that becomes increasingly difficult (in the absence of 
some "tuning" of the soft scales) when the latter are too far above the TeV scale. In that 
case the stability of the electroweak (EW) scale under these quantum effects, with well above 
TeV-scale spartners may become questionable. A way to see this more quantitatively is to 
examine the stability (or relative variation) of the EW scale under (relative) variations of the 
UV parameters (masses) of the theory, with all experimental constraints taken into account. 
This relative variation is the so-called EW fine-tuning measure introduced long ago [3] and 
is a quantitative measure of supersymmetry as a solution to the hierarchy problem. 

A highly stable EW scale under the aforementioned quantum corrections is indeed prefer- 
able, which indicates that the fine tuning measure should be minimal with respect to (vari- 
ation of) these UV parameters and this fact will enable us to make predictions. The corre- 
sponding region of the parameter space can then be regarded as the most natural, and its 
phenomenological predictions are worth a careful investigation. Adding to this analysis the 
requirement that the results found satisfy the WMAP constraint [1] will enable us to put 
together large distance (dark matter) and short distance (EW/TeV scale) physics, and thus 
to improve the consistency and the predictive power of such study. 

The results of such an investigation, at two-loop leading-log level, are reviewed below, 
in the context of constrained MSSM. As we shall see later, the need for a two-loop level 
calculation instead of a 1-loop one is well motivated. This, together with a careful account of 
the dependence of A on the MSSM parameters, will bring interesting, new results. We then 
go beyond the MSSM framework, by considering the effects on the EW scale fine-tuning and 
Higgs mass from "new physics" that may exist beyond this model and which we parametrize 
using effective operators. For technical details of the results we present see [H [2]. 

One can argue that we are missing a quantitative indicator of what upper level of fine- 
tuning is still acceptable or of what makes a model more natural than another. We shall 
not address the latter issue, but for a given model one can accept that, while a given value 
of A may be a subjective criterion of naturalness, the best A is certainly the one that is 
minimized with respect to the UV parameters of the model, as explained above. This is the 
point of view that we adopt and this has some support from the Bayesian method to data 
fitting [2] in which A, as defined in [H] and used below, is automatically present [HIE]- Indeed, 
in this method 1/A emerges naturally as an effective prior in the analysis of a precisely 
measured observable, in this case mz (the EW scale). In this way the prior and therefore the 
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Bayesian method impose automatically a fine-tmiing "penalty" for those points of larger A, 
which are in this way excluded from the final fit. This shows the physical meaning of A and 
the Bayesian preference for points of small A. We shall then search for points of minimal A. 
To place this discussion on more quantitative grounds consider the MSSM scalar potential 

V = ml\Hi\^ +ml\H2\^ - {ml Hi- H2 + h.c.) 

+ (Ai/2) \Hi\'^ + (A2/2) li^sl^ + A3 + A4 \Hi ■ 



+ 



(A5/2) {Hi ■ H2Y + Ae \Hi\' {Hi ■ H2) + A7 \H2\' {Hi ■ H2) + h.c. (1) 



The couplings Xj and the soft masses receive one- and two-loop corrections that for the MSSM 
are found in [HI [9]. The fine tuning of the EW scale with respect to a set of parameters p 
is 13] 
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A _ max I Apl^^^^g ,„2,^2,s2,„2^^}, Ap _ -^j— (2) 

in a standard notation for the constrained MSSM parameters p. 

As mentioned earlier, one could ask whether this otherwise widely used formula for fine- 
tuning is the most appropriate. Variations of this definition indeed exist: for example another 
possibility is to use the "quadrature" version of A defined as 

P 

It is generally agreed that within a model, a value of A or A' > 100 or so (i.e. fine tuning 
worse than 1 part in 100) is rather unacceptable. In any case, one would prefer to have a 
minimization of A or A' with respect to the above UV parameters p. Using this idea, we 
identify the regions of the CMSSM parameter space of minimal A that can then be regarded 
as the most natural. This idea can be used to make predictions, such as to find the most 
natural value of the SM-like Higgs mass m/j, which is important. The results given in the 
following are all based on definition ([2]), however we checked that the predicted value of nih 
obtained from minimizing instead A' is the same, which we find interesting. 

Let us present the general idea briefiy, developed further in the next section. In the 
CMSSM, at tree level mh<mz, and even in the presence of quantum corrections, nih is often 
below the observed LEP2 bound [lOj (114.4 GeV), however large values for mh (up to w 135 
GeV) can indeed be achieved. One would like to clarify what value for ruh in this range 
{mz, 135 GeV) is the most natural, in the light of the criterion of minimal A mentioned and 
in the absence of the LEP2 bound, that is actually not imposed in the following. By doing 
so, we explore the whole range of quantum corrections to nih, from their vanishing to their 
largest values, without the prejudice of imposing a cut on their size, as LEP2 bound would 
do. We return to this problem later in the text, when we discuss the dependence of A on 
radiatively corrected nih and the minimum of A for the whole parameter space. 
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One of the two minimum conditions of the scalar potential V in MSSM gives f ^ = — m^/A 
where m? = m?(j>, j3) is a combination of soft masses (^,1^2,3) while A = A(p, j3) is the effective 
quartic Higgs coupling: 
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where A345 = A3 + A4 + A5. Let us examine the first minimum condition for V: v'^ = — m^/A. 
The problem one sees immediately is that with v ~ 0(100) GeV, m ~ 0{TeV) but with 
A < 1 it is in general difficult to satisfy this minimum condition. To this purpose, with v fixed 
to the EW scale, one has to keep m? as low as possible, closer to EW (scale)^, by "tuning" 
the loop-corrected coefficients entering in the first eq in (jj]) so that the near-TeV scales 
present there largely cancel together, to leave m ~ 0(100) GeV. Another way to phrase this 
problem (sometimes called "little hierarchy" ) is how to separate the EW and supersymmetry 
breaking scale (m^), while still respecting the minimum condition v'^ = —m?/X, with TeV- 
scale soft masses. This "tension" is only one aspect captured by the fine-tuning measure A, 
and is usually more emphasized compared to a more important one, related to the size of A, 
discussed next. 

Obviously, any increase of A will reduce the aforementioned tension and ultimately will 
reduce A. This increase can be due to quantum corrections to the couplings A,, or to other 
corrections to them from "new physics" . This stresses the importance of quantum corrections 
to A in the MSSM. As for the "new physics" beyond the MSSM, this can be represented by 
new gauge interactions, additional effective couplings due to integrating out some massive 
states, etc. Note that the smallness of A, fixed at the tree level by SM gauge interactions, is 
partly to blame for this fine-tuning problem, in addition to the larger values of soft masses 
(^^1,2,3)5 from negative Susy searches or other effects (like demanding the largest quantum 
correction to the Higgs mass to respect the LEP2 bound). Thus the ("little hierarchy") 
problem discussed here is not only one of mass scales (EW versus Susy breaking scale) but 
also of the smallness of couplings Aj. 

We can summarize this discussion by remarking that there are two competing effects, 
quantum corrections or "new physics" contributions that can increase (or mf 2 3) against 
those that could increase effective A. It is indeed possible that the latter corrections dominate 
in some region of parameter space, to allow low fine-tuning, and we shall identify this region. 
Ultimately, if one keeps increasing the soft masses to very large values, these will dominate 
the fine-tuning, to the extreme case of recovering the SM case and the hierarchy problem. 

There is a second minimum condition of V that can be written in the form 
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This condition induces an implicit dependence of tan/3 on parameters p (via A = X{p,f3)), 
that must be carefuhy taken into account when evaluating A |T]; one can then trade Bq for 
tan P which can then be regarded as an independent parameter. With this information, one 
can show that the fine tuning has a general formula [11] (see also the Appendix of [12j): 
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This is the most general dependence of A on the parameters p of the model and has the 
advantage that once the radiative corrections to the couplings and masses are introduced, 
can easily be evaluated. Using available 2-loop leading log corrections for couplings entering 
the MSSM Higgs potential [HI |^ one ends up with a general expression for fine tuning which 
can be used to evaluate A at two-loop. For technical details on this matter see [1]. 

We would like to stress the strong effects that quantum corrections to couplings Aj have 
in reducing the fine tuning. Note that in the limit these corrections to Aj are turned off, one 
can show [T] that after some calculations A of ([HD reduces to the so-called "master formula" 
of [13] (see also [HUS]), which is 
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The numerical discrepancy between the result given by ([7]) and the more general formula 
([6|) is indeed significant and ([7]) often gives larger A than ([6]). Unfortunately, most works in 
the literature use ([7]), albeit with loop corrections to the effective couplings later included in 
the potential, but not in A; this lead to significant overestimates of the overall fine-tuning 
amount. This problem can be avoided if using SoftSusy [16] when the whole procedure can 
be done numerically and agrees well with ([H]) PP in the 2-loop leading log approximation. 

To illustrate the discrepancy between the two expressions, it can be shown that including 
only the one- loop correction 5 from stop/top Yukawa coupling to our Higgs coupling A2 — >• 
^2(1 + (5) and then using ([6|), one finds 
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Since usually 5 = 0(1), one sees that a factor of 2 reduction in A is easily achieved by 5 
(effect not captured by ([7])) even in this limiting case of only one Yukawa correction. More 
corrections to the couplings Aj are likely to reduce A further. Such effects were not included 
in the previous estimates of fine-tuning, and are likely to reduce the CMSSM overall amount 
of fine-tuning, as we shall see shortly. 

There is another way to illustrate the strong impact of quantum corrections on A. Usually 
A ~ mg ~ exp(m^/m^), where the first step comes from first eq in ([3]), while the second 
from the fact that leading quantum corrections to nih are of the type m| ~ logmg, which 
is then inverted into an exponential. Therefore A depends exponentially on the quantum 
corrections to m/j. Even though 2-loop leading-log corrections to rrih may be small relative 
to the 1-loop ones, when "exponentiated" they can bring a significant impact on A. For this 
reason, it is advisable to include not only 1-loop but also 2-loop corrections to Aj. 

In the light of the exponential behaviour of A wrt m/j, the very existence of a global 
minimum of A situated at the intersection of two such exponential dependences on nrih (see 
later) that follow the two minimum conditions, cannot be stressed enough. The results 
reviewed below always use the general formula eq.(l6|) with 2-loop leading-log corrections 
to Aj. 

2 ruh from minimal fine-tuning and dark matter consistency. 

With the above observations, we review the implications of minimal A for the whole parameter 
space. Associating this with the most natural regions of parameters values, we can see 
what this criterion predicts for the Higgs sector. The numerical results for A include two- 
loop corrections with the theoretical constraints: radiative electroweak breaking (EWSB), 
non-tachyonic sparticles masses (avoiding colour and charge breaking (CCB) vacua), and 
experimental constraints: bounds on superpartner masses, electroweak precision data, h — )■ 
s'y, Bg ^ /x"*" n~ and anomalous magnetic moment (5a^. For further details and for the actual 
experimental values considered see Table 1 in [1]. Let us stress that, unless stated otherwise, 
consistency of nih with the LEP2 bound (114.4 GeV [TD]) and/or consistency with the dark 
matter abundance (of the LSP), are not imposed. Later we discuss separately their impact 
on the value of predicted nih. 

The results on fine tuning presented below were obtained using a combined analytical 
and numerical approach, using 2-loop leading log corrections implemented as in [1]. The 
results obtained were in very good agreement with those found using SoftSusy code, whose 
very long CPU run time was reduced (from 6-years on 30x3GHz), using a specially designed 
Mathematica code. This was one reason that prevented earlier investigations of A of similar 
accuracy. The Mathematica code was used to select the phase space points that respected the 
above constraints and also to identify points of minimal A, which were then used in SoftSusy 
and this reduced the run time to manageable levels. 
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Figure 1: Left figure: Fine tuning vs Higgs mass, in a two- loop analysis, for a wide range of parameters 
IJ,o,mo, Aq, Bo,mi/2 and for 2 < tan/3 < 55. The solid line is the minimum fine tuning with central 
values for strong coupling (a3(mz)) and top mass (to*): (as^rrit) = (0.1176, 173.1 GcV). The dashed 
line corresponds to (a3,mt) = (0.1156, 174.4GeV) and the dotted line to (0.1196, 171. 8GeV), to 
account for la variations of and top mass [T7]. The LEP2 bound of 114.4 GeV is indicated by 
a vertical line. Note the steep (« exponential) increase of A on both sides of the minimum value, 
situated near the LEP2 bound for m^,. Right figure: The plot displays, for any fixed p, max|Ap|, 
p = /Xg, jtiq, Aq, Bg , m^^2i ^fiat contribute to overall A of the left figure. The largest of these all, for 
all p and rrih gives the boundary contour presented in the left figure. 

The results obtained are plotted in Figure [TJ In the left figure, 2-loop A is presented as 
a function of the Higgs mass. Interestingly, there is a minimum of A ~ 8.8 which predicts 
a value of nih which is just above the LEP2 bound, at nih = 114 it 2 GeV. The quoted 
theoretical uncertainty of ±2 GeV is due to higher order perturbative corrections that account 
for differences between the results of SoftSusy and FeynHiggs codes and can be even larger, up 
to 3 GeV. This value of ruh is found with the above experimental and theoretical constraints, 
but without imposing the LEP2 bound. It just turns out that minimal fine-tuning wrt to 
UV parameters indicates a size for the quantum corrections that prefer a total value for rrih 
close to this bound. 

Notice the presence in Figured] of the steep, ~ exponential increase of minimal A on both 
sides of its minimum value situated near the LEP2 bound, effect largely due to quantum 
corrections (note the log OY scale) . This behaviour underlines the importance of the minimal 
value of A. Variations of la around the central values of a^{mz) and of the top mass, give 
the results shown by the dotted and dashed lines: a la increase of a^{mz) gives very similar 
results to a 1 cj decrease of top mass; thus these have opposite effects. Indeed a larger top 
Yukawa helps a radiative EWSB, while a larger strong coupling a^{mz) has an opposite 
effect, at two-loop. The minimum of A is where these effects are balanced, in other words 
- assuming minimal A, the Higgs mass prediction cannot be larger due to QCD quantum 
effects. In other words, QCD does not "like" a larger Higgs mass, unless one is prepared to 
accept the fine tuning cost that comes with it (this is indeed very high for rrih > 126 GeV, 
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Figure 2: Left plot: Two-loop fine tuning vs Higgs mass with the influence of the WMAP bound. 
The blue (darker) points sub-saturate the relic density. The red (lighter) points correspond to a relic 
density within the 3a bounds of D.h^ = 0.1099 ± 0.0062 2]- The 'strips' of points at low Higgs mass 
appear due to taking steps of 0.5 in tan/3 below 10. A denser scan is expected to fill in this region. 
Similarly, more relic density saturating points arc expected to cover the wedge of sub-saturating points 
at ~ 114 GeV and A > 30. The continuous line is that of minimal electroweak A of Fig.[l] without 
the relic density constraint. Right plot: as for left plot, but within la. 

when one already has A > 1000). The conclusion is that, from a 2-loop evaluation of A, one 
finds, rather intriguingly, that minimal A ~ 8.8 favours the value rrih = 114 it 2 GeV. 

In the right plot of Figure [J are shown the individual contributions max | Ap|, p = HqjUIq, 
Aq, Bq, rriy^j to the electroweak fine-tuning A. At low m/j, below the LEP2 bound A^2 
is dominant, while above this bound and at large trih, dominant, with A^2 reaching 

similar values near 120 GeV; this contradicts common claims in the literature that A, ,2 
is actually the largest and dominant part of A, (true only below the LEP2 bound). The 
transition between the two regions is happening at about 114.5 GeV, shown also in the left 
plot. With A^2 largely related to the EW effects while A^g related to QCD effects, one 
sees again the significance of the minimal value of overall A, at the interplay of these effects. 
Again, the LEP2 bound is not imposed at any time. 

In principle A^2 (/ij is the top Yukawa coupling) could be included in the overall definition 
of A of eq.©. However, its contribution to A is always sub-dominant (relative to A^g) 
when assuming the modified definition of Ap [20j, appropriate for measured parameters. 
Under this modified definition one must replace: Ap — )• Ap x cTp/p, where Cp is the la error 
of experimentally measured p, in this case hf. With this modified definition Af^2 does not 
change A. See also figure 2 in [1] where A^2 is actually computed and shown but without 
the modified definition. The largest Ap for all rrih in the right plot of Figure [1] generates the 
lower continuous curve of A in the left plot in Fig [TJ 

Let us now put together the short distance (EW/TeV scale) physics effects that we have 
discussed so far with the large distance physics (dark matter) effects. We therefore analyze 
the dark matter constraints on CMSSM, whose LSP is a good dark matter candidate. To 




Figure 3: Left: The 2010 Atlas (CMS) observed exclusion limit given by the red (back) curve in 
the (?Tio,TOi/2) plane, for tan/3 = 3, = and /i > |18| . In our plots of fine tuning in the plane 
(too, mi/2), shown in figures I¥|5I6[ the corresponding red and black exclusion curves are also displayed 
in similar colours. Right: The 2011 Atlas observed exclusion limit (red curve) In our plots of 

fine tuning in the plane (too,toi/2) shown in Figures HJ O El this 2011 exclusion curve is displayed in 
green, to avoid confusion with the exclusion curve from the left plot. Note these exclusion curves are 
for fixed values of some CMSSM parameters {tan/3, Aq, etc). 

this purpose, one takes the phase space points in Figure [H and evaluates for each of them 
the relic density, Qh"^, using micrOMEGAs [21] and test if it is consistent with WMAP [3]. 
The result is presented in Figure [21 where we show the points consistent with the WMAP 
value as well as those that saturate it within 3a. The plot is very similar for a la saturation 
of WMAP value, with only minor differences. Requiring minimal A and consistency with 
WMAP leads to the predictions: 

mh = 114.7 ± 2 GeV, A = 15.0, (consistent with WMAP bound). 

ruh = 116.0 ±2 GeV, A = 19.1, (saturating the WMAP within la). 

mh = 115.9 ± 2 GeV, A = 17.8, (saturating the WMAP within 3a). (9) 

We checked that using a different definition for A, such as A' of ([S]) does not change these 
results for nih, since min(A') is found at similar values for nih and its plot as a function of ruh 
with its dark matter constraints is indeed very similar (not shown here). The main difference 
is an overall shift of the fine-tuning plots of Figure [2] towards higher values of fine tuning, by 
a factor between 1.5 and 2. The minimum of A' indicates however the same value of m/j. 

To conclude, minimizing the fine-tuning together with the constraints from precision 
electroweak data, the bounds on Susy masses and the requirement of the observed dark 
matter abundance lead to a prediction for m/j, without imposing the LEP2 bound, that is 
marginally above this bound. This is an interesting result, and represents our prediction [1] 
for the CMSSM lightest Higgs mass based on assuming A as a quantitative test of Susy as a 
solution to the hierarchy problem. 
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Figure 4: The points in the plane {mQ,mi/2) that have A < 100 and are consistent (blue) with the 
WMAP constraint (Sct deviation) or saturate it (red) within 3cr. These points are the same as in the 
left plot of Figl21 The points in lighter (darker) red/blue have rrih below (above) the LEP2 bound 
for m/j. The black, red and green curves correspond to exclusion curves from CMS, Atlas 2010 curve 
PE] and Atlas 2011 curve [12] respectively, see figure [3] for details. One sees that experimental data 
already test and rule out some points of low fine tuning A < 100. 

While minimal values of A are preferable, one can nevertheless ask what the bounds on 
the parameter space are, for a A beyond which supersymmetry is considered to fail to solve 
the hierarchy problem. It is in general agreed that this happens for A > 100 or so. Therefore 
the bound A< 100 together with the dark matter consistency (3fT upper limit), generates the 
following upper limits on the CMSSM parameters and m/j: 

rrih < 121 GeV, fi < 657.2 GeV, mo < 3.2 TeV 
127.6 GeV < mi/2 < 599 GeV, -1.76 TeV < Aq < 2.26 TeV (10) 

These values can be re-evaluated for a different upper value of A. Note the value of that 
corresponds to EW fine-tuning of 1 part in 100. Given the exponential increase of A with 
rrih, one sees that at mh = 126 GeV one already has a very large, unacceptable fine-tuning 
A = 1000. The mass limits in Tabled] scale approximately as a/A~, so they may be adapted 
depending on how much tuning one is willing to accept. 
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Table 1: Upper mass limits on superpartners (GeV) for which A < 100 (no dark matter constraint). 
If any of these states have masses larger than those shown, will require a fine-tuning worse than 1%. 

It is interesting to compare the results in eq. (|10p for the parameter space of CMSSM 
with A< 100 and dark matter consistency, with the recent exclusion limits of the Atlas/CMS 
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Figure 5: As for the plot in Figure H but with: A < 1000 (top left), A < 200 (top right), A < 50 
(bottom left plot) and A < 20 (bottom right plot). Note that A < 1000 implies m/i < 126 GeV, see 
Figure [T] For A < 10000, the plot is very similar to that for A < 1000 with the only difference that 
the areas near mo = 5 TeV and any mi/2 extend up to 10 TeV. See also the plot in Fig.H]for A < 100. 
The points in dark red/blue satisfy the LEP2 lower bound for rrih, while those in light red/blue do 
not. This bound is imposed at 111 GeV to account for the theoretical error of 2-3 GeV mentioned 
earlier, see text. The black, red and green curves correspond to exclusion curves from CMS, Atlas 
2010 curve [TH] and Atlas 2011 curve [12] respectively, see figure [3] for details. 

experiments \18\ 119] . shown in Figure [3j As it can be seen from Figures HI O [6l LHC 
experiments already probe points in the parameter space that have fine-tuning A < 100. Note 
however that the experimental exclusion curves are for very specific values of the parameters 
(tan/3, Aq, ....) which are not summed over, and therefore points below these curves that have 
different values for these parameters are not tested or ruled out. 

The situation changes for different values of A that one is prepared to accept, see the 
plots in Figures [5] together with the Atlas/CMS exclusion limits. In particular the case with 
A < 20 selects points with the smallest fine tuning. It should be stressed that the presence 
of fewer points (smaller, shrinking area) in the "moduli" space {mQ,mi/2) does not imply 
anything special for the overall fine-tuning amount, like a larger A (disproved by Figures [5]). 
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Figure 6: Left: Two-loop fine-tuning versus Higgs mass for tlie scan over CMSSM parameters with, 
no constraint on the Higgs mass. This is the same plot as in figure [T] but with different colour 
encoding. The dark green, purple, red and black coloured regions have a dark matter density within 
nh'^ = 0.1099 ± 3 X 0.0062 while the lighter coloured versions of these regions lie below this bound 
il/i^ < 0.1285 (3cr). The colours and associated numbers refer to different LSP structures: Regions 
1, 3, 4 and 5 have an LSP which is mostly bino-like. In region 2, the LSP has a significant higgsino 
component, about 10%. Right: Regions of fine tuning A < 100, summed over tan/3 and Aq with 
nih > 111 GeV (to account for the 2-3 GcV theoretical error). Same colour encoding as in the left 
plot. The CDMS-II bound [36] is applied and reduces the area of purple points. The black, red and 
green curves correspond to exclusion curves from CMS, Atlas 2010 curve [H] and Atlas 2011 curve 
[in] respectively, see figure [3] for details. 

It can for example be related to the quality of the scan of the parameter space. Apart from 
this, it indicates the most likely values of the parameter (moduli) space that are preferred by 
the low energy physics, that a fundamental theory (like a string model) should explain or fix 
dynamically. For related interpretations of the recent LHC results see [231 [23] . 

It is important to note that most of the region with A < 100 of the parameter space is 
being tested by the combined LHC (at 7 TeV) and CDMS-II and Xenon experiments [36], 
as it was discussed in |22j together with the possible detection signals. Here we only make 
some comments on the impact of LHC and dark matter experiments on the various regions 
of the fine-tuning plots shown so far. In figure [6] colour encoding corresponds to the structure 
of the LSP and shows where points of different (A,m/i) are located in the plane (mo,mi/2)- 
Red points (1) have a low gluino mass, the LSP is mostly bino with a cross section off nuclei 
too small to be probed by the next generation of direct dark matter searches. Green points 
(3) have lighter squarks, mostly bino LSP and together with the red points are reachable by 
LHC run I (right plot). Black points (5) have gluino and squarks near the TeV scale and 
an almost pure bino LSP. Purple points (2) have a significant higgsino component (10%), a 
heavy gluino, (~ 900 GeV or larger) and TeV-scale squarks. Although they are not within the 
reach of LHC run I, they are sensitive to direct dark matter searches; the CDMS II bound is 



11 



applied in their case while XenonlOO (not applied) can reduce their area further [36]. Also, as 
seen from the left plot the purple region includes points of very low fine-tuning. In conclusion 
dark matter and LHC run I searches are rather complementary in covering the entire plane 
(?TT-o,mx/2) in the TeV region. To cover all the parameter space, beyond A < 100 will require 
running at the full 14 TeV CM energy. 

3 The impact of "new physics" beyond MSSM on A and m/j. 

Following this analysis, some natural questions emerge. How do the above predictions for 
m/i and A change under the presence of "new physics" that may exist at some high scale 
(few TeV or so), beyond the MSSM Higgs sector? Recall that the MSSM Higgs potential 
is the most minimal construction allowed by supersymmetry, but new physics in this sector 
can exist: for example extra Higgs doublets, gauge singlets, additional massive U(l)' bosons, 
which can all affect the MSSM Higgs sector, its quartic effective coupling and its predictions. 
Another question is the following: assume for a moment that a Higgs particle is not found 
at the minimal A as discussed above, is it then possible to have a larger m/j, of say 121 GeV 
but with A ~ 0(10) instead of the corresponding value found above of A = 100? 

One answer is that the model considered is too constrained, and this may lead to a large 
A. Indeed, it is known that gaugino universality condition considered here, if relaxed, reduces 
the value of A [25]; another possibility is to relax the universal Higgs mass which can also 
reduce A. Finally, another possibility is that "new physics" missed by the CMSSM Higgs 
sector, can increase the effective quartic Higgs coupling which leads to a reduction of the 
fine-tuning amount. The ways to achieve this are numerous but also model dependent. 

A simple possibility is to consider the case of MSSM with a low supersymmetry breaking 
scale {\/y ~ few TeV ) in the hidden sector [26\ I27j . In this case, when integrating out the 
auxiliary field of the goldstino superfield that is coupled to the MSSM, one generates apart 
from the usual soft terms, additional quartic terms, without introducing new parameters in 
the visible sector of the model. Indeed, the scalar potential contains a term [27] 

Vd j^\mj\hi\^ + ml\h2\'^ + Bohi.h2f (11) 

in addition to the MSSM potential, in the standard notation. This term can be significant 
for low / while in the more familiar cases of high scale Susy breaking is strongly suppressed. 
Such term can help to increase the SM-like Higgs mass at the tree level and reduce A since 

can be reduced compared to its MSSM counterpart in eq.®, due to the extra term in the 
denominator. For more details of this class of models see |27j and references therein. 
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Other possibilities to increase m/j and reduce A can exist, due to "new physics". To 
perform a model independent analysis, below we consider an effective theory approach, that 
parametrizes in a general way the "new physics" that may exist beyond the MSSM Higgs 
sector. This is done by using a series of effective operators, whose consequences for A and 
m/i are explored below. Let us mention that introducing such operators comes at the cost of 
having new parameters in the theory, beyond those of the MSSM. 

The power of such an effective approach resides in its organizing principle that relies on 
an expansion in inverse powers of the scale M of the "new physics" that generated these 
operators. To this purpose we consider the effective operators of dimensions d=5 and d=6 
that can exist in the Higgs sector. The reason of considering both classes is twofold: first, 
the latter can be present as a leading contribution even in the absence of the former. For 
example integrating out a massive U(l)' generates d=6 operators, with no d=5 ones. The 
second reason is that, if generated by the same physics, at large tan/3 the d=5 operators 
receive additional suppression and then become of similar order of magnitude to the d=6 
operators. In other words, the convergence of the expansion (in 1/M) for large tan /3 requires 
one include both classes of operators, when originating from same massive state that was 
integrated out. 

3.1 The case of d=5 operators and their corrections to mh. 

We start with the d=5 operators in the MSSM Higgs sector. There are two of them: 

= Yl j "^'^ ^'^^^^ {H2.Hif + h.c., X'h{S)/M = Cio + Cii S, (Cio, Cii ~ 
jT^ = ^ J d^e ^ a(5, S^) H2 e"^^] [c(5, S^) e^^ + /i.e.} (13) 

where S = OOrriQ is the spurion field and we assume ttiq <^ M so that the expansion is 
convergent. Here a, 6, c are all dimensionless functions of the spurion that have the form 
a{S, S'^) = aQ + aiS + a\ S"^ + S\ and similar for b and c. 

£1 can be generated by integrating out a massive gauge singlet or SU{2) triplet. Let us 
consider the former case. Indeed, in the MSSM with a massive gauge singlet, with an F-term 
MS^ + E H1.H2, when integrating out S via the eqs of motion, one generates £1. This result 
is similar to considering a generalised NMSSM with a supersymmetric mass term MS^, in 
the decoupling limit of the singlet. Therefore our results for fine tuning are relevant for the 
generalised NMSSM model, in this limit. 

Regarding C2, this can be generated in various ways (see Appendix A, B in [28]) but per- 
haps the simplest way is via an additional pair of massive Higgs doublets of mass of order M. 
However, £2 is actually "redundant", since it can be removed by general spurion-dependent 
field redefinitions, up to soft terms renormalisation, /i term redefinition and 0(1/M^) correc- 
tions [28] . Therefore C2 plays no role in the following. Note that £,2 also includes a particular 
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Figure 7: Left figure (a): the fine tuning A as a function of rrih (GeV) at one- loop; A of MSSM is 
plotted in light blue (mt — 174.4 GeV) with an orange edge (shift induced for rrit = 171.8 GeV) and 
extends up to rrih ~ 114 GeV from where it grows exponentially; A of MSSM with d=5 operators 
with (2/xoCio) = 0.07, (2moCii) = is plotted in dark blue {nit = 174.4 GeV) and with a red edge 
(if rrit = 171.8). Right figure (b): similar to figure (a) but with Cio ~ 0, (2moCii) = —0.1. Non-zero 
or larger Cio,ii (dark blue and red areas) shift the plots to higher rrih, for fixed A. Similar behaviour 
is present for simultaneous non-zero value for both coefficients Ci- 

d = 5 operator of the form J cP9HiOH2 which can be re-written as a Kahler term (for details 
see |33] and the Appendix in the first work in [2]). 

Due to Ci the scalar potential acquires additional terms which bring a correction to the 
effective quartic coupling: A — A + (2CioAio) sin 2/3-1- (— l/2)moCii sin^ 2/3. So A may increase 
and as a result we expect that fine-tuning can be reduced. The correction to the MSSM 
lightest Higgs mass is then 



+ 
+ 



(2Cio^*o)t'^ sin 2/3 



1 ± 



w 



+ 



(-2Ciimo)?^^ 



5mlH, with 5mlH = 0{l/M^) 



It 



m,|) cos^ 2/3 



w 



(14) 



and 



("ii)MSSM 



2 Clo Ato •t- 
sin2/3 



2moCii 



V 



5m\ = 0(1/M2 



(15) 



for the pseudoscalar Higgs. 5m\ ^ and 5m\ are 0{1/M'^) due to d = 5 and, if present, 
also d = Q operators. The upper (lower) signs correspond to h (H), and w is given by 
w = {m\ + m\)'^ — 4 m\ cos^ 2/3. With this result one can show that the mass rrih can be 
increased above the LEP2 bound, for low tan /3, also with the help of quantum corrections. 
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Let us now see the amount of fine-tuning as a function of the Higgs mass, in the presence 
of d=5 operator^. Unhke the MSSM case the analysis of A is done at one-loop only, for a 
sample of points in parameter space with: 1.5 < tan/3 < 10, 50GeV < mQ,mi2 < 1 TeV, 
130 GeV < /Uo < 1 TeV, -10 < At < 10 and 171.8 < mt < 174.4 GeV, consistent with 
the signs for CiOjCii chosen so as to reduce the fine tuning. One-loop analytical results for 
A were obtained (for details see [12]), and the corresponding numerical results are shown 
in Figure [71 Note that in these figures the structure apparent at small A and large m/j is 
a scanning artifact, with the under-dense wedge shaped regions to be filled in with a more 
dense parameter sample. For a fixed value of A and relative to the CMSSM case, one may 
see a systematic shift towards higher nifi. This is due to an increase of the effective quartic 
Higgs coupling (which also increases m/j relative to its MSSM value) and, as a result, also 
decreases A. The overall result is that the minimum amount of fine-tuning A in the presence 
of d = 5 effective operators is reduced relative to the MSSM case, so values of rrih as large as 
130 GeV can still have a low A < 10. The reduction in the fine tuning at low tan/3 (< 10) 
relative to the MSSM case is actually much more marked than that shown, given that in the 
MSSM, in this limit and nifi above the LEP2 bound, A actually increases. 

Let us examine the scale of new physics needed for this reduction in fine tuning. One has 

M^^ = ^^^0^,o, (16) 
uo 0.0/ 

With fiQ between the EW scale and 1 TeV, this shows that large values of M are allowed: 
M w 6(9) TeV for /xq = 200(300) GeV, respectively. Also, in this case A < 10 for 114 < 
iT^h < 130 GeV. To relax these values one can use that an increase (decrease) of (2 /io Cio) by 
0.01 increases (decreases) rrih by 2 to 4 GeV for the same A. 

Let us conclude with a remark on fine-tuning in the presence of a massive gauge singlet 
that was assumed to generate the above d = 5 operator, and a comparison to the NMSSM 
case. The above case is less fine tuned than the NMSSM because of the nature of the super- 
symmetric contribution of the gauge singlet to effective quartic higgs coupling. In our case 
A receives a contribution A ~ sin 2/3 (see eqs. ()14p ). while in the NMSSM the corresponding 
contribution of the singlet is further suppressed, being proportional to A ~ sin^ 2/3. As a 
result the NMSSM is more fine-tuned (for a similar oc 2Au^). This difference can be 
further traced back to the absence in NMSSM of a bilinear F-term MS^, that could increase 
the mass of the singlet significantly (and generated the d=5 operator). Correspondingly, the 
increase of A and thus the reduction of the fine-tuning in NMSSM is not as significant as in 
the case discussed above (for a fixed m/i). 

^ There is a subtle point: notice that the RG eqs are not affected at one loop by the presence of a massive 
gauge singlet that was assumed to generate d in the first instance. Therefore this analysis is consistent. 
The situation is more complicated in the presence of SU(2) triplet; however, the Higgs mass increase and the 
reduction of fine-tuning due to it remain true, even though the exact value of A may be different. 
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3.2 The case of d=6 operators and their correction to m/j. 

We can extend the previous discussion to also include the corrections to nih from effective 
operators of dimension d = 6 that can exist beyond the MSSM Higgs sector. Computing 
the analytical corrections to 111^ from individual effective operators is indeed possible. It 
is however difficult to evaluate the fine tuning A in this case, since the operators can be 
generated in various ways by (unknown) states charged under the SM group, that affect the 
RG flow of couplings in the model and thus the fine tuning. This is unlike the case of Ci in 
the d = 5 case which was assumed to be generated by a massive gauge singlet. Nevertheless, 
it is obvious that an increase of nih by the effective operators reduces the fine tuning A. For 
these reasons, in the following we restrict ourselves to computing the corrections to nih- 
The list of d = 6 effective operators which are polynomial in fields is |12^ [29] 

01 = ^ y d^O Zi {Hi H,)\ = ^ j d^eZ^{H\ e^^ Hi)H2.Hi + h.c. 

02 = ^ y d^O Z2 {Hi H2f, Oq = ^I d^e Zq {H^ e^^ H2) H2.H1 + h.c. 

03 = ^ y d^e Z-i{H\ ifi) {Hi H^), Ot = ^ j dH Z-j Tr W {H^ Hi) + h.c. 

04 = ^ y d'^eZ^ {H2. Hi) {H2. i/l)^ Os = ^j d^'eZs {H2 Hi)^ + h.c. (17) 

where W"' = {-I / 4)1D^ e'"^ D'^ is the chiral field strength of SU{2)l or U{1)y vector 
superfields Vw and Vy respectively. Also Vi^2 = ^?^J(<7"/2) + {^1/2)Vy with the upper sign 
for Vi. Finally, the wavefunction coefficients are spurion dependent and have the structure 

{l/M'^)Zi{S,S^) = aio + aiiS + a*imQSS^ + ai2mlSS\ ~ l/M^. (18) 

where S = mo99. 01.2,3 can be generated in the MSSM with an additional, massive ^7(1)' 
gauge boson or SU{2) triplets, when these are integrated out [30]. O4 can be generated 
by a massive gauge singlet or SU{2) triplet, while Os.e can be generated by a combination 
of SU{2) doublets and massive gauge singlet. O7 is essentially a threshold correction to 
the gauge coupling, with a moduli field replaced by the Higgs, difficult to generate in a 
renormalisable theory with additional massive states. Finally, Og exists only in the non-Susy 
case, but is generated when removing the d = 5 derivative operator C2 by field redefinitions 
[28] . therefore we keep it. There are also operators which involve derivatives (see later). 

Let us see the implications of these operators for the corrections to the Higgs spectrum. 
They bring 0{1/M'^) corrections denoted 5mf^ ^^drin?^ in eq. (fH|) . (fT5]) and their exact expres- 
sions can be found in [21 [31]. However, for most purposes, an expansion of these in 1/ tan (3 
is accurate enough. At large tan/3, d = 6 operators bring corrections comparable to those 
of d = 5 operators. The relative tan/3 enhancement of 0{1/M'^) corrections compensates 
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Figure 8: The corrections Snih due to Oi (left plot) and O3 (right plot) operators, in function of 
the 2-loop (leading-log) rrih of the CMSSM, with M — 8 TeV. Only the Susy part of these operators 
is considered, labelled by aio (0130), respectively. Light blue: CMSSM points with relic density 
nh'^ > 0.1285; on top, in dark blue are CMSSM points with fi/i^ < 0.0913 (3cr deviation); on top, in 
red, are CMSSM points that saturate WMAP bound within 3cr (WMAP value nh"^ = 0.1099 ±0.0062). 
The total, corrected value of Higgs mass is then rrih + Smt- The CMSSM points below the lower 
continuous line have A< 100, those between the two continuous lines have (100< A<200) while those 
above the upper continuous line have A > 200. Therefore 5mh for A < 100 is rather small, of few 
GeV. O2 {O4) gives results remarkably close to those of Oi (Os), respectively, and are not shown. 




Figure 9: As for Figure [51 but for the supersymmetric part of Os (left) and (right) operators. 
Note the minus sign in front of 5mh showing that for the chosen (positive) sign of parameters (ajo) a 
decrease of rrih is actually obtained. 

for their extra suppression relative to 0{1/M) operators. However, in some models only 
d = 6 operators may be present, depending on the details of the "new physics" generating 
the effective operators. If niA is kept fixed, one finds in ©(l/M^) order [2] (see also [3T]): 
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a22"io + (aso + a4o)/io + 2a6i ?wo/^o -a20™| - (2 Cio ^o)^ ("t-a - "^1)"^ 

+ t;^cot/3 (m^-m|)"^(^4m^ ( (2a2i+a3i+a4i+2a8i)"io^o + (2a50+a6o) /^o + a62?n-o) 

- (2a6o-3a7o) m|- (2a6o + a7o) ?7t-|) + 8(m^+m|)(/iomoCioCii) ^^^/I'TT-i-^^l)^ 
+ 0(l/tan2/3) (19) 

The full value of is that obtained by using eq. (fT9]) in (fn|) where (w-^) mssm is replaced by 
the 2-loop leading log valu^. So the effective operators correction is regarded as a classical 
correction ("perturbation") added to the 2-loop leading-log CMSSM value, and crossed-terms 
that involve products of loop-corrections and effective operators coefficients are neglected, 
being of higher order. For an easier reading of the effects of the operators, one can introduce 
the correction dm^ 



Sruh 



1/2 , 1 5m 



"^^l2-loop,CMSSM + - "^'^l2-loop,CMSSM = 0^;^ + 0(1/M^) (20) 

"'^l2-loop,CMSSM 



with dmi as in (fT9|) . The total Higgs mass value is then 5mh + "2,ft|2-ioop,CMSSM- Fur- 
ther, it is preferable to search for possible increases of m| by supersymmetric rather than 
super symmetry-breaking effects of the effective operators, because the latter are less under 
control in the effective approaclJ§. In any case, the non-Susy part of the effective operators 
has an impact on 5mh that is comparable to that of the supersymmetric part considered here. 
So in the following we concentrate only on Susy corrections, induced by the coefficients ajo 
with j to label the corresponding operator Oj. 

The results obtained are illustrated in Figures [8l [9] where the correction 6mh is shown 
as a function of the CMSSM value of m^ (2-loop, leading-log). The results are obtained in 
the the following way. We consider all the phase space points displayed in Figure [2] that 
respect all experimental and dark matter constraints (except the LEP2 bound on mh that 
is not imposed). On this "background" we consider the perturbation due to the effective 
operators of dimension d=6 and evaluate 5mh as outlined above and shown in these figures, 
as a function of the 2-loop leading-log value of mh in CMSSM. 

From these figures one can conclude that the CMSSM points with lowest fine-tuning 
A < 100 corresponding to m^ < 121 GeV, have a rather small 5mh, of few GeV. For M = 8 
TeV dmh is up to 4 GeV (6 GeV for A < 200), and this decreases (increases) by 1 GeV 
for a 1 TeV increase (decrease) of AI. However, the value of M is usually restricted to 
be in the region of 8 TeV or higher, from p parameter constraint^ [32] and our expansion 



^In the numerical evaluations and figures below we used the exact expression of [2], see also [31| . 
■^Also, one would prefer a supersymmetric solution to the fine-tuning problem associated with increasing 
the MSSM Higgs mass well above the LEP2 bound. 

''This bound applies to a combination of operators and can in principle be reduced for individual operators. 
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parameter m?aij was always taken < 1/4, where m is any scale of the model {mQ, mi/2 



n). This is considered conservative enough for a convergent perturbative expansion. Note 
that points which were below the LEP2 bound by the correction 5m/j mentioned, become 
now phenomenologically viable. The special point of CMSSM of minimal A = 17.8 that 
saturates the relic density within 3a and with irih = 115.9 it 2 GeV, could receive a correction 
6mh ~ 4 GeV, so that m/j increases to rrih + 5mh = 119.9 it 2 GeV. Correspondingly, its 
earlier associated A is likely to decrease by a factor proportional to the square of the ratio 
of the final to the initial value of ruh- 

Given the relatively small size of their correction 5mh one can say that the particular 
CMSSM points with A < 100 and their predictions are rather stable against the presence 
of supersymmetric "new physics" at the scale M = 8 TeV. This finding can be explained 
by the fact that these points generically have a light /j. and also light rni/2, and thus the 
supersymmetric corrections 6mfi, (oc fJ-'^ajo, etc) are rather suppressed. The corrections Snih 
can increase or decrease if one also includes effects of Susy breaking associated to Oi, but 
these bring additional model dependence and extra parameters. 

Let us now discuss about the CMSSM points in Figures [H [9] with fine tuning A > 200, 
situated above the upper continuous line. They can bring an increase of nih which can be 
significant, of 10-30 GeV. For example there are points which for nih near 100 GeV can receive 
corrections of order 15-20 GeV, to reach and comply with the LEP2 bound. Interestingly, for 
Ci^2 the Higgs mass increase is such that total m/j remains close to ~122 GeV. Points that 
are largely fine-tuned and have a value for uih significantly below the LEP2 bound, are often 
receiving the largest corrections Snih- This opens the possibility that the phase space of the 
CMSSM be increased and points which were otherwise ruled out on grounds of extreme fine 
tuning and/or LEP2 bound, can be "recovered" and become phenomenologically viable. 

Finally, notice that we kept all operators Oj independent of each other. By doing so, one 
can single out the individual contributions of each operator (labelled by Oj*), which helps in 
model building, since not all operators are present in a specific model. Also it is unlikely that 
"new physics" will bring up in the leading order, simultaneously, all these operators. What 
does all this mean for EW scale fine-tuning? The value of A for those points initially strongly 
fine-tuned can decrease by a factor equal to the square of the ratio of the Higgs mass after 
and before adding the correction 5mh and this effect can be significant, as seen earlier for Ci 
{d = 5 case). A similar effect is expected for the case of d=6 operators. So one expects a 
change of A 



One cannot obtain a more exact evaluation of A in our model-independent approach, i.e. 
in the absence of the details of the new physics (quantum numbers of massive states) that 
generated the effective operators in the first instance. For further discussions on this see [2]. 



771' 
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(21) 
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3.2.1 Effective operators: removing redundant operators. 

The reader may notice that the discussion for the d=6 operators ignored a class of operators 
that could be present (and that are similar to Lq, of the d=5 case). These operators involve 
extra derivatives and are 

09 = ^ y d^Q ^9 V' V2 Oi2 = ^l (fe Zi2 Hi V" W^^) H2 

010 = ^ y d^O ZiQ Hi V' H2 = ^ j d'^O Z13 h\ e^^ W^^^ V" Hi 
On = ^ y d^e Zn if| e^^ V" W^^^ Hi Ou = ^ J d^O hI e^^ wj^^ V° H2 
Oi5 = ^ J d^OTie^ We-^ D^e^ Wa e"^) (22) 

Here VaHi = e~^^Dae^^Hi and is the field strength of Vi. O15 does not depend on 
Higgs explicitly, but could affect its scalar potential. To be general, in the above operators 
one should include spurion (5) dependence under any Vq to account for supersymmetry 
breaking effects associated to them. Such operators are easily generated when integrating out 
massive states in 4D models. They can also be generated at one-loop by compactification, 
after integrating out towers of Kaluza-Klein states, in the presence of localised (Yukawa) 
interactions [3l] (for example Og, Oio) or due to bulk (gauge) interactions [35] (Ois). 

Such operators are however redundant in the sense that they can be removed by non-linear 
field redefinitions [21 [28]. To see how this works, consider for example operator Og, without 
the gauge field dependence, and in the presence of an otherwise standard Lagrangian, with 
an arbitrary superpotential W: 

. ^ (1 . o/M^) * + ,{J.^S Wl.^. xl . ..cj . 0(VM=) (23, 

Here ^^n^ comes from Cg with the replacement Hi — )• $. Further, one replaces ^^n^ _). 

(-l/l6)I^^t/)2^ rpj^-g £ 

can be unfolded into a "standard" Lagrangian without extra 

derivatives (see [33], also Appendix B in the first paper in [2]). To see how this works, 

consider a change of basis to $1.2: ^ = si^i + 82^2 and (1/m) D^^^ = ri^i + r2^2 where 

si,2, fi,2 form an unitary matrix, so that the eigenvalue problem is not changed; m is a 

very small, non-zero mass scale of the theory that can be taken to zero at the end of the 

— 2 

calculation. Since <^ and D <1> are not independent, such transformation must be accompanied 
by a Lagrangian constraint, which must vanish in the limit Af — )■ 00. This constraint has the 
form: 

6C = y"d20[(l/m)D^(si^>i + S2^'2)^-(n^>i + r2$2)](mV(4M))$3 (24) 
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where is a chiral superfield which plays the role of a Lagrange multiplier, so that the total, 
equivalent Lagrangian in the new basis is CJ = £ + 5£. One then brings CJ to a diagonal 
form to find the result: 



+ 



T^[«>(l>i,2);x] -Meal's 



1 



M3 



(25) 



where $ = <I>2 — $1 and where we took the limit m — t- 0. This Lagrangian is that of a 
second order theory, only polynomial in superfields, and classically equivalent to the initial 
one (j23|) . Given the signs in the D-term, two massive superghosts (^2,3) are present, of mass 
near the effective cut-off, 0(M). Their presence is just an effect of truncating the operators 
series expansion to l/M^ termjf|. Also observe that the x field was "spectator" throughout 
this analysis and did not affect it; in fact the x-dependence can be replaced by an arbitrary 
polynomial function. Further, since the superghost degrees of freedom are massive, one can 
integrate them out, by the equations of motion. After a careful calculation and consistent 
Taylor expansion, the result is (see appendix in [2]) 



1 



^\^^-^'^'\^r^xW\'^^\x\ + x\\^\ I d^dW[^i;x]+h.c.}+0 



r 1 



(26) 



where the derivatives of W are taken wrt its first argument. This Lagrangian contains only 
polynomial interactions (renormalisable or not) and standard kinetic terms, and is equivalent 
to the original one, eq. ([25|) . but has no extra-derivative terms. 

This result agrees with that obtained by using the equations of motion in the derivative 
term in (p3]) . but this is not true in general, as we argue below. Indeed, as shown in the 
appendix of [2], the presence of a derivative ternj§ $n$/M in the superpotential of an 
otherwise arbitrary Lagrangian leads, via the method shown above, to a result that is similar 
in order 0{1/M) to that obtained via the eqs of motion. However, in the 0{1/M^) the result 
found by using ordinary eqs of motion is different and actually wrong. The discrepancy is due 
to the fact that Euler-Lagrange eqs are changed in the higher order theory, and this should 
be taken into account when using eqs of motion to eliminate the higher-derivative operators. 
We find the method presented above more elegant and transparent. 

Returning to the CMSSM with an additional Og, the second term in ()26p together with the 
standard MSSM superpotential (~ fiHi.H2) can bring only wavefunction renormalization of 
the Higgs kinetic term or other effective operators polynomial in fields. When Susy breaking is 
included, such terms also bring soft terms and /x-term redefinition. A similar strategy applies 



® Indeed, in a renormalisable or ghost-free theory with a massive state, when integrating out this state (by 
the eqs of motion) one finds in the low-energy effective action an infinite series of derivative terms, suppressed 
by powers of M. The "truncated" theory will have a finite number of derivatives, and, as seen above, this will 
bring superghosts (or just ghosts in the non-Susy case), i.e. fields with negative kinetic term. This happens 
even though the theory with the infinite series is ghost-free, since the original theory was so. 
Such term is part of £2 discussed in the case of d=5 operators since 
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to the rest of the operators hsted in (j22p . Therefore all these operators are "redundant" 
in the sense that they can be eliminated, up to redefinitions of the fields, soft masses and 
//-term |28j . For this reason they were not considered in the phenomenological studies of the 
previous sections. 

4 Conclusions 

A quantitative test of supersymmetry as a solution to the hierarchy problem is the fine-tuning 
A of the electroweak scale with respect to variations of the UV parameters of the model, after 
including the quantum corrections and experimental and theoretical constraints. While the 
largest value of A for which Susy is still a solution to this problem is somewhat subjective, 
it is very natural to ask that A be minimal in a given model. This point of view is also 
supported in part by Bayesian approaches to data fits, in which A is automatically included 
in the effective prior expression as an extra 1/A factor, so that the method brings naturally 
a fine-tuning penalty for points of large A. This underlines the physical meaning of A in 
general and the need of minimizing the overall fine-tuning as done in this work, to select the 
points of physical relevance. 

We applied this idea of minimizing A to the CMSSM at two-loop leading log, and inves- 
tigated its results for the value of the radiatively corrected value of the lightest Higgs mass, 
rrih- Remarkably, although A depends ~ exponentially on nih, there does exist a minimum 
at a very acceptable value (A = 8.8) corresponding to a mass m/j = 114 ± 2 GeV. The very 
existence of such a minimum situated at the intersection of two exponential dependences on 
both sides of this value for ruh, and induced by quantum corrections, cannot be stressed 
enough. The exponential growth of A for nih above this value is largely due to QCD quan- 
tum effects which overcome the "good" Yukawa loop effects needed to induce radiative EW 
symmetry breaking. Thus QCD "does not like" a larger nih unless one is prepared to pay the 
associated (high) fine-tuning cost. Imposing consistency with the WMAP bound, the above 
value changes mildly to A = 15 leading to nih = 114.7 it 2 GeV, while saturating this bound 
within 3(7 leads to A = 17.8 with nih = 115.9 =b 2 GeV (the quoted theoretical uncertainty 
(±2 GeV) can actually be larger, up to 3 GeV). 

It is indeed remarkable that constraints from short distance physics (EW precision data) 
and from large distance physics (dark matter) can be consistent with each other so accurately, 
and together can help one to make physical predictions. We also checked that using a different 



definition for A, such as A' = yX^jA? does not change this result for nih, since minA' is 
found at similar values for nih and its plot as a function of nih is very similar. Finally, let 
us note that points with A < 100 are currently being tested by the LHC 7 TeV run and by 
dark matter experiments (CDMS, Xenon, etc) and their impact was briefly discussed. 

We further discussed how the results for A and mh change under the addition of "new 
physics" beyond the MSSM Higgs sector and also whether one could have a larger m-h with 
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a low A. This is relevant for MSSM since for irih > 121 GeV, A is already A > 100, 
and becomes 1000 for nih = 126 GeV! The "new physics" could reduce A to acceptable 
values, even at larger m/j. It could be represented by additional Higgs doublets, massive 
gauge singlets or massive U(l)', etc, and its effects can be generally parametrized by series 
of effective operators, suppressed by the scale of massive states that generated them. We 
considered both d=5 and d=6 operators. In the presence alone of the former, one can 
achieve a low A < 10 even for nih as large as ~ 130 GeV. Such a d = 5 effective operator 
can be generated by integrating out a massive gauge singlet beyond MSSM, whose effect is 
to increase the quartic higgs coupling and reduce the fine tuning for similar ruh- This case is 
nothing but the decoupling limit of the generalised NMSSM which contains a supersymmetric 
mass term for the gauge singlet (in ordinary NMSSM, the fine tuning reduction is smaller). 
Other ways to generate such operator can however exist. 

Regarding the operators of dimension d=6, we considered their individual, supersymmet- 
ric corrections to the CMSSM 2-loop leading log value of m/j, for an expansion parameter 
m/M < 1/4 with rh any scale of the model (mo, /n-1/2, ^, ^^)- Their effects were treated as a 
perturbation of the CMSSM "background" points that respected all experimental constraints. 
It was shown that CMSSM points with A < 100 having ruh < 121 GeV (A < 200, ruh < 122 
GeV) receive supersymmetric corrections of up to 4 (6) GeV from individual operators, re- 
spectively. Therefore, the points in the CMSSM of lowest fine-tuning receive rather modest 
(supersymmetric) corrections from individual operators, and are therefore rather stable un- 
der "new physics" at a scale considered here of 8 TeV. Applied to the point of minimal fine 
tuning (114 ± 2GeV), the correction mentioned of 4(6) GeV from individual d=6 operators 
brings ruh up to a value of 118(120) ±2 GeV. Including a 3a consistency with WMAP changes 
this result mildly up to 119.9(121.9) it 2 GeV. Finally, an increase (decrease) of the scale M 
by 1 TeV brings a decrease (increase) of the correction by about 1 GeV. A next step in this 
study is to impose dark matter constraints on the CMSSM with d=6 operators, since by 
supersymmetry these extend the neutralino sector too. The obtained bounds on M can then 
be used to re-evaluate the quoted correction to m/j in the presence of these operators. 
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